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Abstract 

We study the process of relaxation back to thermal equihbrium in (1 + l)-dimensional 
conformal field theory at finite temperature. When the size of the system is much larger 
than the inverse temperature, perturbations decay exponentially with time. On the other 
hand, when the inverse temperature is large, the relaxation is oscillatory with character- 
istic period set by the size of the system. We then analyse the intermediate regime in 
two specific models, namely free fermions, and a strongly coupled large k conformal field 
theory which is dual to string theory on (2 + l)-dimensional anti-de Sitter spacetime. In 
the latter case, there is a sharp transition between the two regimes in the k = cxd limit, 
which is a manifestation of the gravitational Hawking-Page phase transition. In particu- 
lar, we establish a direct connection between quasinormal and normal modes of the gravity 
system, and the decaying and oscillating behaviour of the conformal field theory. 
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1 Introduction 



An important problem in finite temperature field theory and statistical mechanics is to study 
the response of a system in thermal equilibrium to a generic perturbation. In particular, one is 
typically interested in understanding the process of relaxation back to thermal equilibrium. For 
small perturbations, this is well described by linear response theory, and the time evolution of 
the relaxation is determined by the retarded correlation function of the perturbation [H I2l E] • 
Generally, in the presence of interactions this problem is addressed within finite temperature 
perturbation theory. A special role is played by scale invariant theories, where the zero tem- 
perature 2-point functions are uniquely determined (up to a normalisation) by scale invariance. 
However, finite temperature introduces a new scale and consequently the conformal Ward iden- 
tities no longer determine the Green's functions completely. However, significant progress can 
be made if the conformal field theory (CFT) has a dual formulation in terms of gravity (or 
string theory) on anti-de Sitter space (AdS) spacetime. Indeed, the AdS/CFT jH [HI El Ej 
duality predicts that the retarded CFT correlation functions are in 1-1 correspondence with 
Green's functions on anti-de Sitter space with appropriate boundary conditions [HlinillllllllllE]- 
Furthermore, the poles of the retarded CFT correlators are given by the quasinormal modes 
in AdS jHl E] • This correspondence was verified explicitly in the high temperature regime of 
the 2-dimensional CFT, dual to supergravity on AdSs [TT]. Other calculations of quasinormal 
modes in anti-de Sitter backgrounds have appeared in [13]. 

The purpose of the present paper is to analyse finite volume effects for a quantum conformal 
field theory living in one space-like dimension (closed to form a circle of length L). In two 
limiting cases, when the dimensionless parameter LT is infinite or zero, conformal invariance 
completely determines the correlation function, independent of the details of the theory in 
question. The behaviour of the retarded correlation functions in these two cases is qualitatively 
different. In the first case, the perturbation decays exponentially with characteristic time 
proportional to the inverse temperature. In the second case, we have oscillatory behaviour 
with characteristic period determined by L. Our main purpose here is to study the relaxation 
process in the intermediate regime when LT changes from zero to infinity. 

After describing the qualitative behaviour based on general arguments, we analyse quan- 
titatively the case of a non-interacting theory, and contrast its behaviour with the strongly 
coupled large k CFT dual to supergravity on AdSs (see also [Ellin])- (The parameter k plays 
the role of N in the usual terminology of large N CFT dual to AdS gravity.) In the latter case, 
we have an explicit expression for the 2-point function at finite temperature and finite volume. 
We can thus analyse the linear response of the CFT. In the limit k = oo, there is a sharp 
transition between a regime of exponential decay and oscillation. This is a manifestation of the 
gravitational Hawking-Page phase transition between the BTZ black hole and thermal AdS. 
Moreover, we then establish a direct connection between the behaviour of the linear response, 
and the behaviour of the corresponding bulk perturbations. In particular, we show that the 
regime of exponential decay is governed by the quasinormal modes of the BTZ black hole, while 
the regime of oscillation is governed by the normal modes of thermal AdS. Thus, the behaviour 
of the bulk AdS perturbation is mirrored precisely by the behaviour of the linear response of 
the CFT. Finally, we speculate on the expected behaviour of the strongly coupled CFT at finite 
values of k. 
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2 Linear Response Theory 



We consider a system initially in thermal equilibrium, and apply a small perturbation. The 
main goal of linear response theory is to study the change in the expectation value of an operator 
0{x,t) as a result of this perturbation. The total Hamiltonian of the system takes the form 
H'{t) = H + Hcxt{t), where H is the unperturbed Hamiltonian, and Hext{t) couples an external 
field to the system, with the assumption that i^ext = 0, for t < to. The change in the ensemble 
average of O is given by 



6{0{x,t)) = i[ dt' Tr{p [H,xt{t'),0{x,t)]} , (1) 
where p is the unperturbed thermal density matrix. If we take the perturbation to be 

H,xt{t) = J dx J{x,t)0{x,t) , (2) 
where J{x, t) is the external source, then (Q) takes the form 

dt' / dx' J{x', t')D^{x, t; x', t') , (3) 

-oo J 

where 

t; x', t') = -ie{t - t')Tr{p [0{x, t), 0{x', t')]} (4) 

is the retarded propagator. 

For the particular case of an instantaneous perturbation of the form 

J{x,t) = 5{t)e"''' , (5) 

one finds 

^e-''^'D\uJ, k) . (6) 
-oo 27r 

In the following, it will be useful to recall the Lehmann representation of the retarded propagator 

j-oo Iti uj — uj + le 

where the spectral function p{uj,k) is the Fourier transform of the commutator p{x,t; x' ,t') = 
{[0{x,t),0{x',t')]). 

From general arguments, we know that D^{uj, k) is an analytic function in the upper half 
tu-plane. If the energy levels of the system are discrete and the ensemble sum converges, then 
we conclude that D^{uj, k) is a meromorphic function of uo in the lower half plane, with simple 
poles on the real axis. Correspondingly, as can be seen from the linear response will show 
oscillatory behaviour. Indeed, as shown in PHl-CH]! the energy spectrum for a system with 
finite entropy is discrete, and will in general exhibit the phenomenon of Poincare recurrence. If 
the energy levels are continuous, the properties of D^{uj, k) can be more complicated. At zero 
temperature, the retarded Green's function will generally have poles and cuts on the real axis. 
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corresponding to stable states and multi-particle states, respectively. Furthermore, there can 
also be poles in the lower half cu-plane corresponding to resonances. The distance of the poles 
from the real line then determines the decay time of such a resonance. At finite temperature, 
D^{lj, k) is a meromorphic function in the lower half plane. Eqn. ^ then shows that the 
characteristic times for the thermalisation of an instantaneous perturbation is determined by 
the imaginary part of the poles of D^{uj, k) in the complex cj-plane. For a generic interacting 
theory, the location of these poles will depend non-trivially on the coupling constants, and 
concrete results are known only in specific limits, where perturbation theory is applicable. 



3 Conformal Field Theory 

To see how these general results are realized in (1 + 1) -dimensional conformal field theory, 
we first calculate the correlation functions on a torus with periods L and (3 = T~^. The 
correlation function in real time, t, is then obtained by the analytic continuation r = it, where 
r is imaginary time. When the size of one direction on the torus is taken to infinity, the 
torus becomes a cylinder. Using Cardy's result j2Q|, the correlation function on the cylinder 
is obtained from the correlator on the plane by the conformal mapping^ w = [l/^ir) In^, w = 
(Z/27r) Inz, where w and z are complex coordinates on the cylinder and the plane, respectively. 
The 2-point function on the cylinder is then given by 



[smh['K / l)[w — w')J^"[smh(7r//)(iL' — w')Y'^ 

where (/i, h) are the conformal weights of O. 

Suppose now that the temperature is finite, while the size L is taken to infinity; then 
/ = ± yu), where /i is the chemical potential and w = a + it. After analytic continuation 
r = it, we have w = a — t and w = a + t. We then obtain the real time correlation function 



D±{t,a) = (O(tT«e,a)O(0,0)) 



[sinh TTTnia - t ± ie)]2'^[sinh 7rTi((T + t t ie)] 



2h 



(9) 



where = 1// and Tl = 1/1 are the effective right and left temperatures. Clearly, the correlation 
function decays exponentially for large time t] for example, D±(t,0) ~ ^-2nT{h+h)t ^ when Tl = 
Tr = T. 

For {h, h} G |N, the Fourier transform of the commutator, p = (D^ — D^), can be evaluated 
by the method of residues leading to 



p{uj, k) oc 



2'kTrJ V 2nTL 



(10) 



where p± = ^{u ^ k). From (fTUj). we see that p{uj,k) has an infinite set of simple poles on 
either side of the real line. The poles of the Green's function D^{uj, k) are then obtained by 

^Note that in 2-dimcnsional CFT the left and right sectors are independent. In particular, they can be 
defined on cylinders with different radii. 
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restricting this set to the lower half plane. These are given by ^T] 



uj = k — Ani Tiln + h) , 

u = —k — Arci Tfi{n + h) , n G N . 



'Ill 



On the other hand, if the size L is kept finite and the inverse temperature (3 = is taken 
to infinity then / = L. The direction of the imaginary axis can be chosen in the a direction, 
w = T + ia. The analytic continuation to real time now results mw = i{t + a) and w = i{t — a) . 
Thus, we find that (jH)) leads to the real time correlation function 



(7r/2L)2(^+^) 



[(sin l{tTie + a)]2'^[(sin ^(t ^ k - a)] 



2h 



(12) 



This clearly shows the expected oscillatory behaviour in real time. In the following, we specialise 
to the subset of operators with h = h & N /2. To isolate the pole structure of the retarded 
Green's function, we again consider the Fourier transform of the commutator 



p{uj,k) 



dt 

-oo Jo 



da e 



iujt—ikcr 



(13) 



with k = ^m. We first consider the case with h 



h = 1/2. Then p takes the simple form 



p(o;,fc)ocE^(^-^ 

nGZ \ 

Thus, we find the poles 



n 



al^^ I 
u m 

2n ' 



2n 



\m\ 



(14) 



2tt. . An. 

— \m\ + —{n + h) , 

1j Lj 

271. . An. 

\m\ [n + h) 

L L 



(15) 



with m G Z, and n G N. The spectral function for higher values of h takes a form similar to 
(fT^ . with additional prefactors designed such that the pole structure is given by (fl^j) . 

The considerations so far are quite general, and rely only on the conformal properties of the 
operators driving the perturbation. The behaviour of the correlation functions and (fT2|) in 
the two limiting cases is thus universal, and in agreement with general expectations. However, 
the behaviour in the intermediate regime, when both L and T are finite, cannot be derived 
entirely from conformal symmetry. In order to understand this regime, we first consider the 
2-point function of free fermions on the torus |2I]. This takes the form. 



(^(W7)^(0)), 



Ou (f I ^|) d^e, (o I ) 

(O I 9, (l I ^L 



(16) 



where w = r + ia and v characterises the boundary conditions for tp{w). For finite temperature 
boundary conditions, we have z/ = 3,4. Using the properties of 6'-functions, it is then easy to 
see that (jl6j) is invariant under shifts w w + (3 and w ^ w + iL. 
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The real time correlation function is obtained from by the substitution w = i{t + a) . 
The resulting real time correlator is thus a periodic function of t with period L. Zeros of 
the theta function 9i{w / P\iL / P) are located at [21j w = m/3 + inL, where m,n are integers. 
Therefore, for real time t, the correlation function ()16|1 is a sequence of singular peaks located 
at (t + cr) = nL. In the limit L/ (3 oo, the correlation function (fTT)|) approaches (actually 
the left-moving part of with h = 1/2) which exponentially decays with time. In the opposite 
limit, when L/(3 — > 0, it approaches the oscillating function (fT^ . This is in agreement with 
our discussion using Cardy's arguments. It is important to observe that the exponential decay 
of the function in the limit of large L is consistent with periodicity of the real time correlation 
function (fTHIl with period L. In order to see this, we note that the correction to the leading 
behaviour is governed by e-'^L/f^e'^i^+'^)//^ _ Thus, the limit of large L is meaningful only for times 
much smaller than L. For such t, the exponentially decaying function is a good description of 
the correlation function. However, for t approaching L the correction terms become important, 
and the time-periodicity of the correlation function (|TB|l is restored. 



4 Strong Coupling Regime 

At strong coupling, conformal symmetry is not enough to determine the linear response at 
finite volume. We are therefore unable to obtain exact results in the general case. However, 
progress can be made in some special cases. Here, we consider the supersymmetric conformal 
field theory dual to string theory on AdSs. This theory describes the low energy excitations of 
a large number of Dl- and D5-branes ^ [7[ El- It can be interpreted as a gas of strings that 
wind around a circle of length L with target space T'^. The individual strings can be simply or 
multiply wound such that the total winding number is k = |, where c » 1 is the central charge. 
The parameter k plays the role of N in the usual terminology of large N CFT dual to AdS 
gravity. In order to obtain information about the correlation functions in the strong coupling 
regime of this theory, we can appeal to the AdS/CFT correspondence. According to this 
duality, each supergravity perturbation $(m,,s) of mass m and spin s propagating on AdSs has a 
corresponding operator 0(^f^ ^ in the dual conformal field theory. This operator is characterised 
by conformal weights {h, h), with h + h = A, h — h = s, and A is determined in terms of the 
mass of the field. The CFT correlators are then determined in terms of the corresponding bulk 
Green's functions. In addition, there is a correspondence between (quasi)normal modes in the 
gravitational background and the poles of the retarded Green's functions in the conformal field 
theory [S]. 

According to the original prescription jUj , each AdS space which asymptotically approaches 
the given two-dimensional manifold should contribute to the calculation, and one thus has to 
sum over all such spaces. In the case of interest, the two-manifold is a torus (t, cr), where /? 
and L are the respective periods. There is an SL{2, Z) family of AdSs spaces which approach 
the torus asymptotically [T31 122] • For the purpose of understanding the Hawking-Page phase 
transition, it suffices to consider the BTZ black hole and thermal AdS space, corresponding 
to anti-de Sitter space filled with thermal radiation. Both spaces can be represented [23 as a 
quotient of three-dimensional hyperbolic space H^, with line element 

ds^ = ^7 (dzdz + dy"^) , y> . (17) 

y2 

The BTZ black hole (for simplicity we consider only the non-rotating BTZ black hole) has 
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inverse temperature /? = 27r/r+, where r_|_ is the horizon radius. Also, z ~ e?"^^/!^ and y ~ e27r(T//3^ 
where w = a + it. Thus, the orbifold identification is given by 

^^e^-^/^z, y^e'^'^'fy. (18) 

For the thermal AdS at the same temperature we have z ~ e^™/^ and y ~ e^'^^''^, where 
w = r + icr. In this case, the identification is 

In both cases, the boundary of the three-dimensional space is a torus with periods L and 13. In 
particular, the two configurations (thermal AdS and the BTZ black hole) are T-dual to each 
other, and are obtained by the interchange of the coordinates r ^ a on the torus. 

The correlation function of the dual operators then contains the sum of two contributions 
as 

{0{w, w)0{w\ w;'))Torus = e-'^^^ {O 0')btz + e'^^^^^ {O 0')ms , (20) 
where 

5'btz = ' 'S'Ads = ^^^^ ' (21) 

are the Euclidean actions of the BTZ black hole and thermal AdSa, respectively [15j. We 
see that these contributions are dual to each other under interchange of (3 and L. Using the 
AdS/CFT correspondence, the correlation functions on the torus were computed in |2^l25] for 
the case of a scalar field {h = h). For the BTZ background, the result takes the form |23 

(O(u..^.)O(0,0))btz = E ,...i,|.(.. + „L)]4nh||(^ + nL)]^> ' 



where w = a + ir. Note that ()22|) takes the form of the strip expression (jSj) summed over images 
to make it doubly periodic. On the supergravity side, the justification to sum over images comes 
simply from the fact that the correlator solves a Green's function equation. From the CFT 
point of view, this result is non-trivial. Indeed, for a generic CFT, summing over images does 
not produce the correct finite volume correlator. For example, the free field correlator (fT^ does 
not have this form. 

Expression (j^UI) is the result for the correlation function in the strong coupling regime. 
Although each term in the sum ()20|) is not modular invariant, the sum over the full SL{2, Z) 
family does have this property fHl 122] ■ Depending on the ratio L/(3, one of the two terms 
in (j2Uj) dominates |15.]. For high temperature {L/(3 is large) the BTZ is dominating, while 
at low temperature {L/j3 is small) the thermal AdS is dominant. The transition between the 
two regimes occurs at /? = L. In terms of the gravitational physics, this corresponds to the 
Hawking-Page phase transition [211 • This is a sharp transition in the limit k = cxd, which is the 
case when the supergravity description is valid. In this limit, the BTZ black hole is the sole 
contribution for L > (3, while thermal AdS is the only term which contributes for L < (3. 

The two terms in (j2m) have drastically different behaviour as functions of real time. After 
the analytic continuation r = zt, the BTZ contribution (f^^ produces the correlator 



(0(t, a)O(0, 0))btz = E + + ^ + ' (23) 
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which exponentially decays with time. Furthermore, in the limit L/j3^ oo, only the n = 
term in ()23p contributes, and we thus recover the universal behaviour 

It is clear that the Fourier transform of the spectral function for ()23|) is again given by 
fjlO|) . with k restricted to the discrete values ^m. Thus, the poles of are given by (fTT|) with 
k = ^m. At first sight, this appears to be in contradiction with the general behaviour discussed 
in section 2. Indeed, for finite volume we expect a discrete energy spectrum and consequently 
that the poles of should lie on the real axis. The resolution of this puzzle lies in the peculiar 
properties of the CFT under consideration [7]. For L/f5 > 1, the typical configuration consists 
of a relatively small number of multiply wound strings so that the effective volume relevant 
for the energy gap is L^s — kL oo. Of course, this explanation immediately raises another 
puzzle: If the effective volume of the theory is infinite, how come that the correlation function 
is periodic in a with period L. The reason for this lies in the fact that the operator O does 
not distinguish between simply wound and multiply wound strings. Consequently, O still sees 
a finite volume. 

For h ^ h, which corresponds to fields of non-zero spin in AdSs, the finite temperature 
Green's functions have not been worked out in the literature. However, it was shown in (IJ 
that there is a 1-1 correspondence between the poles of the retarded finite temperature Green's 
function and the quasinormal modes for fields of spin s = h — h'm the BTZ background. Quasi- 
normal modes are solutions to the wave equations which are purely ingoing at the horizon and 
subject to the boundary condition that the current vanishes at infinity^. Only a discrete set 
of modes satisfying these boundary conditions is possible, and the frequencies are shown in 
[TT] to be identical with (lllj). In ^T], the correspondence was shown for L/ (3 = oo. However, 
the above discussion shows that this is in fact valid for L/f] > 1. Indeed, one finds that the 
bulk perturbation decays via these quasinormal modes in precisely the same way as the linear 
response of the conformal field theory given by ©. 

The question then arises as to the behaviour for L/f3 < 1. On the gravity side, ()2Up 
implies that the thermal AdS Green's function gives the relevant contribution to the real time 
correlation function. The result for the thermal AdS is obtained from (|22p under the interchange 
T ^ a and P ^ L. Hence, we have 

m, 0)0(0, 0))..s = E I,., . ^ ^ ^ , _ ^ ^ . (24) 

which is periodic in t with period L. Clearly, it represents a periodic sequence of singular peaks 
at t±(T = nL. In the limit of infinite (3/L, only the n = term contributes and (|24j) approaches 
the expression (jl2p . The Hawking-Page transition is thus a transition between oscillatory 
behaviour at low temperature and exponentially decaying behaviour at high temperature. From 
the CFT point of view, this behaviour is explained by the fact that in the low temperature 
phase the generic configuration is given by simply wound strings so that the effective volume 
is finite and consequently the energy spectrum is discrete in agreement with general argument. 

For h G |N, the Fourier transform of the spectral function for the finite temperature 
AdS Green's function (j2D) is again given by (HD), with simple poles given by (fT^ . On the 
AdS side, the set of frequencies ffl^j) is identical to the set of normalisable modes these 

^Originally, the quasinormal modes were required to satisfy Dirichlet conditions at asymptotic infinity '10'. 
However, as shown in this leads to problems for h < 1 which can be resolved by requiring the vanishing of 
the current. 
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modes are regular at the origin of AdSs and are normalisable at infinity. As a result, the 
oscillating behaviour of the bulk perturbation is mirrored by the oscillating behaviour of the 
linear response. In j^ZI; the normalisable modes were obtained for a scalar field h = h. However, 
as shown in [2H], higher spin s equations of motion in AdSa can be reduced to that of massive 
scalar fields for any s. We thus expect that the correspondence between normalisable modes 
and the poles of D^{uj, k) in the low temperature phase is, in fact, valid for any s = h — h. 

5 Discussion 

As pointed out in JEj-jn!, the energy spectrum for a system with finite entropy S is discrete, 
with level spacing of the order . As a consequence, the system will necessarily show Poincare 
recurrence. These recurrences will occur on a timescale of the order t ~ . Indeed, as we 
saw in section 2, for a system with discrete spectrum, the retarded propagator has simple 
poles on the real axis. Therefore, the linear response will exhibit oscillatory behaviour. In 
general, however, this oscillatory behaviour will be quite complicated. Typically, one expects 
the correlation function for a system with finite entropy to be a quasiperiodic function, with 
incommensurate frequencies. Poincare recurrence ensures that the evolution is unitary with no 
loss of information. This is exactly the behaviour seen in the free fermion correlation function 
()16|): at finite L it is oscilating with period L. 

Based on these remarks, it is important to understand the behaviour observed in the strongly 
coupled CFT. In the strict k = oo limit, the sole contribution to the correlator is given by the 
BTZ black hole. The decay of this correlator is due to the fact that the effective volume 
Leff ~ kL is infinite, and thus the spectrum is continuous. This explains why we found complex 
poles corresponding to quasinormal modes of the black hole. However, for finite k the correlator 
must ultimately exhibit Poincare recurrence. 

In the present paper, we considered the inclusion of the contribution from thermal AdS. 
We found that this gives a periodic contribution, whose frequencies are the same as the normal 
modes of AdS. While the addition of this contribution does prevent the decay of the correlator 
at late times ^3] , it is not sufficient to produce the Poincare recurrences at finite k. One could 
consider the inclusion of the remaining members of the SL[2, Z) family of solutions with torus 
boundary. These contributions are necessary in order to ensure modular invariance. However, 
they will be parametrically of the same order as the thermal AdS contribution. The total 
correlator will still include the decaying BTZ part with complex frequencies. It seems that in 
order to see the discreteness of the energy spectrum on the CFT side, one will need to include 
finite k corrections to the correlator. In this way we expect that the k = oo (BTZ) contribution 
will be dressed by 1 /k corrections, so that the correlator at finite k will no longer be a decaying 
function of time. One can see how this may happen by recalling the case of free fermions. The 
correlation function is decaying when L = oo and is periodic at finite L. A somewhat 
similar behaviour is expected in the strongly coupled case. As a result, the Poncare recurrence 
time will become finite and set by Lgff. This problem certainly warrants further investigation. 
In effect, there is no fundamental difference between the free and interacting case. The free 
system is periodic in time, while the interacting system should exhibit Poincare recurrences as a 
quasiperiodic function. In both cases, the evolution will be unitary with no loss of information, 
as expected for a system in finite volume. 
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We should also stress that the conformal field theory dual to supegravity on AdSs is very 
special. For a generic CFT, we do not expect a phase transition even at strong coupling. It 
would be interesting therefore to find more examples of interacting CFT's in various dimensions, 
where explicit non-perturbative results can be obtained. In this respect, the duality between 
the 0{N) sigma model in 3-dimensions and fields of even spin in AdS4 might be of interest 
|29j . Note also that while the explicit computation of finite temperature Green's functions in 
AdSd, d > 3, is generally not possible, the quasinormal modes can nevertheless be obtained 
numerically piU l3Uj . In this way, one can obtain quantitative, non-perturbative, information 
about thermal Green's functions also in higher dimensions. For d > 3, the high temperature 
phase is an AdS-Schwarzschild black hole whereas thermal AdS is dominant at low temperature. 
This Hawking-Page transition is then related to linear response theory in the dual CFT via 
quasinormal and normal modes in the two backgrounds. 
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